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Abstract. In this work we study a scenario with a warped 5D smooth braneworld with 4D 
Minkowski geometry builded from bulk scalar matter non-minimally coupled to gravity with 
^f) • an additional Gauss-Bonnet term. We present exact solutions for the full braneworld con- 

figuration in contrast to previous results where only approximate solutions were constructed 
due to the highly non-linear character of the relevant differential equations. These solutions 
^ ' allow us to study the necessary conditions for the finiteness of the 4D Planck mass and the 

■ 4D effective cosmological constant. This fact enables us to perform a more rigorous analysis 

of 4D gravity localization compared to approximate approaches. We also analyze the local- 
ization properties of scalar, vector and tensor fluctuation modes for the constructed field 
configurations. We show that for the considered backgrounds, only the massless tensor mode, 
i.e. the 4D graviton, is localized on the brane, while the vector and scalar modes are not 
confined to the brane. 
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1 Introduction 

Braneworlds models [1-9] are an interesting alternative to the standard Kaluza-Klein 
compactification. In order to protect the standard 4D physics of some unobserved effects, in 
the Kaluza-Klein paradigm extra dimensions are compactified to a tiny size (for a 4D observer 
that is equivalent to the very hight energy regime). The braneworld scenario has another 
approach, in this alternative the physics of our 4D world is compatible with the existence of 
infinite extra dimensions. A procedure to do that is by considering that the standard model 
(SM) fields are trapped on a 4D hypersurface, called a 3-brane (our universe). In contrast with 
the ordinary matter, gravity and exotic matter can reside in the whole higher-dimensional 
manifold (bulk). Although the gravitational field can propagate through all dimensions, one of 
the first requirements to have realistic braneworlds models is to recover standard 4D gravity 
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on the brane. For example, if the induced geometry on the brane is flat we need to localized 
a 4D graviton on our brane. 

The braneworlds models is classified with respect to their width in thin and thick branes. 
In a scenario with a thin braneworld [7-9] the curvature scalars are singular at the location 
of the branes, that is direct consequence of the null width of the branes. In spite of that 
complication, the called Israel-Lanczos junction conditions [10-12] make harmless such sin- 
gularities when we study the localization of gravity and matter and in calculation of several 
effective parameters of the model i.e. 4D Planck mass, etc. 

The fact of considering that the ordinary matter fields are completely confined on a 
region with null width is just an approximation. The phenomenology of the Standard Model 
are well known at electro-weak energy scale mEw^ therefore, although the standard matter 
can not move freely through the extra dimensions, it is possible that they can access to distant 
r < iTi^^ ~ 10~^^ m without contradictions with the Standard Model [5]. In other words, 
there are more realistic alternatives to thin brane configurations they are known, generically, 
as thick braneworlds or domain walls. When the brane has a non-trivial width it is not 
necessary to use the Israel-Lanczos junction conditions, and this might be a mathematical 
advantage whether the action of the model is complicated (too many matter fields, non- 
minimal coupling, higher order curvature terms, etc). 

Thick branes might be constructed by using different procedures. One of the more 
evident example is by replacing the delta functions in the action by a non-singular source 
functions. Another simple method is by using self-interacting scalar fields coupled minimally 
to gravity [13-21]. 

Within the frame of thin (Randall-Sundrum) braneworlds, bulk scalar fields have been 
investigated, for instance, in [22-26], where the Randall-Sundrum scenario has been modified 
by considering self-interacting scalar fields non-minimally coupled to gravity [22]. Further- 
more, in [23] and [24] the resulting action has been modified by adding a Gauss-Bonnet term. 
The influence of higher curvature terms in scalar-field-generated thick brane models has been 
studied in Refs. [27-37]. 

In dimensions higher than four (in 5D, for instance) the standard Einstein-Hilbert action 
may be supplemented with higher order curvature corrections without generating, in the 
equations of motion, terms containing three or higher order derivatives of the metric with 
respect to the space-time coordinates [38, 39]. A particular combination which satisfied that 
requirement is the well-known Gauss-Bonnet (GB) invariant: 

^GB = R'^^'^^'Rabcd - 4.R^''Rab + R^ (1.1) 

where A,B,C,D = 0,1,2,3,5. While in 4D the GB invariant is a topological term which 
does not contribute to the classical equations of motion^, in dimensions higher than four 
the GB combination leads to a ghost-free theory and it appears in different higher dimen- 
sional contexts. For instance, in string theory the GB-term appears in the first string tension 
correction to the (tree-level) effective action [41-47]. In the framework of thick braneworld 
scenarios generated by a self-interacting minimally coupled scalar field, the GB invariant 
has been studied in connection with the localization properties of the various modes of the 
geometry in Ref. [27] (see also, for instance, [36]). The non-minimal interaction between a 
bulk (self-interacting) scalar field and gravity has been considered, for instance, in Ref. [48], 
where the perturbative stability of a thick brane was explored. 

^In scenarios with AdS4 the 4D Gauss-Bonnet term produces non-trivial results in the conserved quantities 
of the theory [40] . 
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It seems quite natural to investigate the effect of a self-interacting scalar field non- 
minimally coupled to gravity [49-51], in thick braneworld scenarios with higher curvature 
terms (a Gauss-Bonnet invariant, for instance). In this work we study a 5D thick braneworld 
modeled by a smooth scalar domain wall non-minimally coupled to gravity with a Gauss- 
Bonnet term on the bulk. Our results are a extension and generalization of previous results 
obtained within the framework of braneworlds with 4D Poincare symmetry generated by 
scalar fields minimally coupled to gravity [52] (by using a similar approach, thick braneworlds 
with de Sitter 3-branes were studied in [53]). In order to accomplish that, we first construct 
exact solutions for the highly non-linear differential equations for the scalar field, a necessary 
step for providing a rigorous study of the consistency and localization properties of the 
metric fluctuations of our braneworld configuration. Moreover, these solutions are essential 
for analyzing the positiveness of the 4D Planck mass (an indispensable property of any viable 
theory). We further analyze the tensor perturbations of the geometry which can be classified 
into scalar vector and tensor sectors according to the 4D Poincare symmetry group. We 
establish that both the scalar and vectorial sectors are not localized on the brane in contrast 
to the 4D massless mode, the 4D graviton, of the tensorial fluctuation modes. We finally 
summarize our results and conclusions at the end of the paper. 

2 The model 

Let us explore a thick braneworld described by the following 5D action (similar set up 
is studied in references [27], [48]): 



where ol > 0, and k ~ 1/M'^, where M - the 5D Planck mass. Besides, is a real scalar 
field and V = V{ip) its self-interaction potential. The quantity L = L{ip) describes the non- 
minimal coupling between the scalar field and the Einstein-Hilbert term and 'T^qb is the 
5D Gauss-Bonnet term (1.1). 

In the scalar-tensor gravity theory (2.1), 15 degrees of freedom qab, plus the scalar 
field propagate the gravitational interaction. Hence, this is not a pure geometrical theory 
of gravity. In particular, the metric coefficients define the geodesic motion of test particles, 
while the scalar field determines locally the strength of gravity by means of the effective 
gravitational coupling oc K/L{(p). 

The way the scalar field ip is coupled to the curvature in the present theory, deserves an 
independent comment. Actually, in this paper, just as a matter of necessary simplicity, we 
have considered explicit coupling of ip to the curvature scalar R, but not to the Gauss-Bonnet 
invariant. Intuition, instead, dictates that the scalar field should couple in a similar fashion 
to R and to TI'q^, since both contribute towards the curvature of space-time. This will be 
the subject of forthcoming work. We want to underline that this 'asymmetric' coupling has 
nothing to do with physical requirements, but just with simplicity of further mathematical 
handling since, as we shall see further, the relevant equations of motion are highly non-linear. 

The Einstein's field equations that come from the action (2.1) take the following form: 




(2.1) 



L Rab + e Qab = K TAB + Va^bL + - qab OL, 



(2.2) 
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where, 

Qab = \ 9 AS ^GB - 2 i? Rab + 4 Rac R^b + ^R^^ Racbd - 2 Racde Rb^^, 

is the Lanczos tensor representing the Gauss-Bonnet corrections to the Einstein's field equa- 
tions, e = 2a' K, □ = g'^^Vc^ D, and the stress-energy tensor tab, corresponding to the 
scalar matter content on the bulk, is defined as usual: 

2 

TAB = dA^dsV - - QAB V{ip). 

The remaining terms in the right-hand-side (RHS) of (2.2) come from the non-minimal 
coupling of the scalar field to the curvature scalar. 

Furthermore, the equation that describes the dynamics of the scalar field (Klein-Gordon 
equation) can be written as follow: 

1 dV 

where, = dL/dip. As it is done, for instance, in Ref. [54], the geometry of our braneworld 
model is described by a warped metric in conformally flat coordinates: 

ds"^ = a'^{w)[r]^^dx^'dx'' - dw"^], (2.4) 

where the variable w is the extra coordinate, and rj^i, is the 4D Minkowski metric. Due to 
the simplicity of our metric ansatz, it is natural to consider that the field (p depends only on 
the fifth coordinate w. By taken into account the last comment and (2.4) the field equations 
(2.2) and (2.3) read 

Sni^ip' 1 



4f-.2 I n/'\ I nil i 01 /2^ 



2Ka? 



i2 



■W [W + W) - {%' + 3^2) L 



(2.5) 



-L^^" -lnL^^' + (i + ^)^''-- {n' -n')q = o, 
/ + m^' - - ^ (sn' + 2n') = o. 

dtp K ' 

Here the tilde denotes derivative with respect to the extra-coordinate, in addition to this 

T-L = —, while q = L ttT-L^. 

a 

As one can straightforwardly check the equations (2.5) are not independent. 
3 Planck masses 

The 4D effective theory can be obtained by integrating the 5D action (2.1) with respect 
to the fifth coordinate w (dimensional reduction). In order to do that let us consider a 
generalization of the 5D line-element (2.4) 

ds^ = a^{w)\gy,y{x)dx^dx'' - dw\ (3.1) 
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where, g^y{x), is an arbitrary 4D metric. The 5D fundamental theory is reduced to a 4D 
Einstein-Hilbert effective action plus the corrections that come from the scalar matter and 
higher curvature terms of the bulk. In this section we will focus in the analysis of the 4D 
Planck mass, therefore we only need to extract the 4D Einstein-Hilbert effective after the 
dimensional reduction of the 5D action^ 

r2 



54 ^ M^i j (^x^/\g^\R^ + • • • , (3.2) 

where the subscript 4 labels quantities computed with respect to 4D metric, g^u{x) and Mpi 
is the effective Planck mass. 

A method to find the low dimensional effective theory and in particular the 4D Einstein- 
Hilbert action consists in interpreting the warp factor as a conformal function as follows [55]: 



9AB H> QAB 

and, rewriting the 5D action in terms of the quantities defined with respect to qab- In order 
to obtain the Einstein-Hilbert part of the 4D effective action, it is only necessary to consider 
the above conformal transformation on the terms, —L(ip) R/2k, and, a'TZ^^, in (2.1). The 
following expressions display these quantities after the conformal transformation. 



R = a-^i^R-8ni3-12{V^fj, (3.3) 

nlj^ = ^ i^n^j^ - i2R(y^f - 24i?Si? + asr^^ (va^b^ - Va^^b^) 

{mf + (W)^ - (VaVbi9)(V^V^i?)1 , (3.4) 



+ 72 



+ 144(VAViji?)(V^i?)(V^i?) + 144(Vi?)2Si?| , 



where, i? = In a. The terms with a tilde are defined with respect to the metric qab- 
The relation between 4D and 5D Planck masses can be written as follow 



L + ^{n^ + 2n') 



dw, 



/oo 
-oo 

/oo 
a^{w)qdw + 8M^e[a']'^^. (3.5) 
-oo 

This quantity should be positive and finite for a consistent theory. Moreover, it should re- 
produce the effective gravitational couplings that we observe in our 4D world, if we wish to 
recover 4D gravity on the brane. Later on we shall see that these conditions are fulfilled for 
a wide class of solutions of our scalar-tensor model. 

4 Background Geometry and Non-minimal Coupling 

Here we shall consider the following coupling function [22] (see also [48]): 

L{ip) = l-^^^. (4.1) 



^Later on we will study the whole 4D effective theory that arises after integrating the 5D action with 
respect to the extra-dimensional coordinate. 
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If the parameter = 0, the scenario is described by a bulk scalar field minimally coupled 
to gravity (it was explored in Ref. [27]). Hence, the parameter ^ characterize models with 
non-minimal coupling between scalar matter and gravity. Positivity of L leads to (p^ being 
bounded from above: 

L > ^ < I ^ ^ Vf • (^-2) 

We shall consider, besides, a regular geometry of the form (2.4), which interpolates between 
two asymptotically AdS^ space-times, which is depicted by the following warp factor [27]: 

a{w) = , °° (4.3) 

where the width of the thick brane is 1/6, and the quantity ao is dimensionless and related 
to the radius of the asymptotic AdS'^ space-time. All of the resulting quadratic curvature 
invariants are regular and asymptotically constant. 

Let us make a change of variable w ^ x = bw where x is dimensionless. By using (4.1) 
the second equation in (2.5) can be rewritten as: 

3 4e62 3x2 



(l + x2)2 al (l + x2)3' 

where, now, the tilde denotes derivative with respect to the dimensionless variable x. It 
is difficult to find the general solution to the above equation, however, several interesting 
(particular) exact solutions can be found. Here we shall split the analysis into three separated 
cases corresponding to: i) a scalar field minimally coupled to gravity and including a 5D 
Gauss-Bonnet invariant (^ = and e ^ 0) [27], ii) a scalar field non-minimally coupled 
to curvature without the Gauss-Bonnet term (^ 7^ and e = 0), and iii) the more general 
situation when both non-linear eff'ects are present, i.e. when ^ 7^ and e 7^ 0. Once a given 
particular (exact) solution for the scalar field ip = ip{x) is found, one hence can write the 
self-interaction potential as a function of the (dimensionless) extra-dimensional coordinate 
x: 

(4.5) 



l + xV 2 

The above is just a rewriting of the first equation in (2.5) for the case of interest in this 
paper. 

Once ip, and V, are given as functions of the dimensionless variable x, the components 
of the stress energy tensor tab = tab{x), are also known functions of the dimensionless 
extra-dimensional variable. 



5 Mirror Symmetry 

An important remark on the particular solutions we are looking for is related to the 
symmetries of the metric coefficients. As mentioned, in the general case when there is a non- 
minimal coupling between the scalar field and the curvature scalar (^^ 7^ =^ L = L{p) 7^ 1), 



- 6 - 



gravity is propagated both by gAB and by if. Hence, one should naively expect that the 
(real) scalar field, ^{x), will respect the same symmetries as the metric functions qab- In the 
remaining part of this section we shall show that, as a matter of fact, while the latter is a 
mandatory requirement for the case when ^ = (in the minimal coupling case), this property 
is not present in the general case in which ^ 7^ 0. 

In the present case the metric coefficients in (2.4) inherit the 'mirror' symmetry 

xi-^-x =^ gAsix) = gAsi-x), (5.1) 

distinctive of the warp factor (4.3). In consequence, the tensors Rab, and Qab, in the left- 
hand-side (LHS) of Einstein's field equations (2.2) 

Rab{x) = Rab{-x), Qab{x) = Qab{-x), (5.2) 

also respect invar iance under (5.1). 

For the minimal coupling case (^ = L = 1), since the Einstein's field equations 

read, 

Rab + ^Qab = ktab, (5.3) 

the invariance (5.2) will entail that, necessarily, TABix) = tab{—x), which means, in turn, 
that 

V{x) = V{-x), ^{x) = ±ip{-x), (5.4) 

where, in the last equality, only one of the two signs, "-|-", or "— ", is to be chosen at once. 
In other words, (f can be either even or odd, under (5.1). 

In the general case ^ / (L = L{{p) 7^ 1), on the contrary, since the coupling function L 
is multiplying the Ricci tensor in the LHS of (2.2), mirror symmetry (5.1), (5.2) is respected 
only if the additional requirement 

L{ip{x)) = L{ip{-x)), (5.5) 

is fulfilled.^ In the case studied in this paper, since L is quadratic in ip (see Eq. (4.1)), the 
latter requirement will entail, again, that, under (5.1), can be either even, or odd. 
In case if does not show any obvious symmetry under (5.1), 

(f{x) / ±(p{-x) L{f{x)) + L{lp{-x)), 

the coupling function L counteracts the symmetry displayed by Rab and, in consequence, the 
LHS of (2.2) - the pure geometrical part of Einstein's equations - is not invariant under (5.1) 
anymore. Hence, if ^ 7^ 0, mathematical consistency of the solutions does not impose any 
mirror symmetry requirement on ^p. This can be only an independent 'ad hoc' requirement. 

Notwithstanding, in the present paper the symmetry requirements (5.1), (5.2), (5.5), are 
be used to select physically suitable solutions among the set of all possible exact solutions we 
will be able to find. Consequently, we shall take into consideration only those (/^-configurations 
which are either even {'p{x) = ip{—x)) or odd {'p{x) = —ip{—x)) under x 1— )• —x. While for 
the minimally coupled case it is legitimate, as shown above, in the case when ^ 7^ 0, it is not 
required by mathematical consistency, but, by symmetric aesthetic, instead. 

^In our analysis we took into account the fact that the operators VaVs, and □ in the RHS of Eq. (2.2), are 
invariant under (5.1). Recall that, in the present case, there is functional dependence on the extra-coordinate 
X only. 
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6 Four-dimensional effective theory and cosmological constant 



In section 3 we obtained the 4D Planck mass in terms of the 5D parameters. As men- 
tioned above, after one performs a dimensional reduction of (2.1) the 5D theory is reduced to 
a 4D Einstein-Hilbert effective action plus the corrections that come from the scalar matter 
and higher curvature terms of the bulk. Thus, the 4D effective theory must be well defined. 
It is not enough to have a finite and positive 4D Planck mass, e.g., the 4D cosmological con- 
stant has to be finite as well, otherwise the model (2.1) will not reproduce the physics of our 
world. Let us study the whole effective theory that arises after the dimensional reduction of 
the fundamental model (2.1). By considering the metric (3.1), the Ricci and Gauss-Bonnet 
term (3.4) take the form 

R= L^R + m' + 12%^), (6.1) 

T^GB = A(^GB + ^2n\n^ + m')) + • • • , (6.2) 

<ijQi 

where the • • • represent terms associated to the 4D Einstein-Hilbert action that were studied 
in section 3. 

Furthermore, R^abc = R^A = 0, thus, the non-trivial components of the Riemann and 
Ricci tensor are the associated to the 4D metric g^i,{x). In order to simplify the analysis let 
us assume that the metric g^i,{x) is a small deviation from ?7^jy, in other words, g^i^i-c) = 
ilfiu + Sgfiui-c), while we can set — t- = + (5(^4 (x), where ip is the fundamental scalar field 
and 5ip4^{x) is its 4D part. The effective action in four dimensions can be written as follows 

Ai 




2 , . ,,,, + X2U4^a) + hT^icBj + Sa, (6.3) 

where the subindex 4 is associated to the metric g^i,{x), S4 is 4D Einstein-Hilbert action 

(3.2) and C/4((^4) is the self-interaction potential associated to the field (p^ = 6ip4{x). On the 
other hand, A4, Ai, A2 and A3 are effective coupling constants and take the following form 

A4 = i y du; (^2a^L{2n' + 3^^) + .^^^ + ^^^y + UeaU'^im' + ^^) 

a^dw, 

a^dw, 

A3 = — / adw. 
3 J 

The low dimensional theory (6.3) is described by a 4D scalar field minimally coupled to gravity 
plus the Gauss-Bonnet term in four dimensions. Furthermore, A4 is the 4D cosmological 
constant. It is not difficult to show that the matter coupling Ai and A2 are finite for the metric 

(4.3) . Similarly to the Planck mass (3.5), the effective cosmological constant A4 depends on 
the fundamental scalar field ip. In the case where \ip\, \(p'\ and \V\ are bounded along the 
extra dimension (these requirements are satisfied for the solutions that we will study in this 
work) the effective cosmological constant is finite. Otherwise, A3 is unbounded. Therefore, 
the effects of the 4D Gauss-Bonnet term are decoupled of the 4D phenomenology. Maybe 
this is not as harmful as it seems, if one takes into account the fact that the contribution of 
the Gauss-Bonnet term to dynamical equations is trivial in four dimensions for almost all 
the spacetimes. 
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7 Exact solutions 

As already mentioned, here we shall split the study into three particular cases i) e 7^ 
and ^ = - minimal coupling, ii) ,^ 7^ and e = - non-minimal coupling, and iii) e 7^ and 
^ 7^ - the general case. We will rely on a number of assumptions, so that, after achieving 
considerable simplification, particular exact solutions can be found. 

7.1 Minimal Coupling Case (e 7^ and ^ = 0) 

If one sets, ^ = 0, then (4.4) can be easily solved. Actually, let us first change to the 
variable 

„2 



r = arctan(x) 4^ x = tan(T), sin r 



2 ^ 



(l + x2)A = i. ^ {l + x^)X' = X. (7.1) 
ax or 

Notice that x G ] — 00, +oo[ r € ] — 7r/2, 7r/2[. In terms of this new variable equation (4.4), 
with ^ = 0, it transforms into the following first order equations: 



3 / 4e62 



which can be integrated in quadratures to give: 



^^{T) = ±^l-E{smT,k) + ^^, (7.3) 

where 

£;(sinT,A;)= [ J 1 - sin^C^C, 
Jo 



k 



'4e62 



is the incomplete elliptic integral of the second kind [56], (p^ are arbitrary integration con- 
stants, and the it signs represent different branches of the solution. 

The particular case when, oq = \/4e6 =^ k = 1 [54], ^^{t) = ±y^3/K sinr + (/?q , or, 
in terms of the dimensionless variable x (see Eq. (7.1)): 



As before, the it signs describe two possible branches of the solution. To fix the constants 
(Pq, one : 
instance, 



, one has to choose appropriate boundary conditions. One simple way is by choosing, for 



Under the latter, specific, boundary conditions at infinity, the solutions (7.4) can be expressed 
as: 

^^{x) = ±\[^(^^-l 

These solutions, however, do not meet the necessary symmetry requirements (5.1), (5.2), 
(5.5). In fact, the scalar field solutions (7.3), (7.4), do not respect invariance under x 1— >■ — x, 
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unless the constants (/Jq are set to zero: ip^ = 0.^ Hence, for instance, the appropriate 
boundary condition at infinity for ip^{x) in (7.4) will be. 



(p^{x — )• oo) = 93, 



0, 



instead of the above written one. 

Summarizing: particular exact solutions of equation (4.4), for the minimally interacting 
case = 0), are: 



^^{x) = ±^l-E 



.k 



^^{x)=±J- 



k^l, 
k = l. 



(7.5) 



These solutions for the scalar field are odd under (5.1), (5.2), (5.5): ^p^{x) = —ip^{—x). 
Furthermore, self-interaction potential V is asymptotically constant and negative (see fig. 
1). This result is consistent with the fact that our geometry (2.4) is asymptotically AdS^. 




w 



Figure 1. The self-interaction potential for the case e 7^ and ^ = 0. In this figure we set k = 1, 
ao = I, b = 1 and n = 1/2. 



7.2 Non-minimal Coupling Case {C ^ and e = 0) 

If we make the same replacements 2; 1— )• r = arctanx (Eq. (7.1)), then the master 
equation (4.4) can be rewritten in the following form: 



Then, without loss of generality, one can make the following assumption: 



cp^ = h{ip) 



1 dh{Lp) 



2 ' 

so that Eq. (7.6) can be written as a first-order differential equation, 

— h 2w = itp- — , 

dip if t^if 



^Recall that, i5(sinr, fc) = — sin r, fc). 



(7.6) 



(7.7) 



(7.8) 
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where 

The fohowing expression for h{f) solves Eq. (7.8): 

h{^) = -^^-^ + C^-'-, (7.9) 
2(1 + uj) t.uj 

where C is an arbitrary constant. 

Here, for the sake of simphcity of mathematical handling, we set C = in Eq. (7.9). 
Then, after substituting back (7.7) into (7.9), one is left with the following first order differ- 
ential equation: 

Eq. (7.10) can be integrated in quadratures: 



* ' ^ r + C„ (7.11) 



2(2^ - k) 



where we have returned to the original parameter ^, and Ci is an integration constant, which, 
in the following calculations we set to zero to meet the imposed symmetry requirements (5.1), 
(5.2), and (5.5). Hence, depending on the interval in ^-parameter, one will obtain different 
particular exact solutions: 

1. < C < k/2, 



or, in terms of the original dimensionless variable (see Eq. (7.1)): 

(/3^(x) = ifQi cos (±/3i arctanx) , (7.13) 
where, for compactness of writing, we have introduced the following constants: 



2(K-20 3^ 
¥'01 = -\/-r7 TT^ Pi-~ 



The bound (4.2) on 93^ leads to: 

, 2(k - 2£) 2 

which, for < ^ < k/2, is always fulfilled. Hence, positivity of L{ip) does not impose 
any additional constraint on the parameter ^. On the other hand, for the above profile 
of the field the self-interaction potential is constant and negative at tt; — )• 00 (see fig. 
2). 
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Figure 2. V corresponding to solution (7.13). In tliis figure we set ^ = 1/10, oq = 1, & = 1 and 

K=l/2. 



2. K, 



or, in terms of x = bw, 

(p^{x) = (/5o2 cosh(ib/32 arctanx), (7.15) 

where. 



In this case, positivity of the coupling function L{ip), imposes the following non- 
algebraic constraint on ^: 



erf ?.LJ1_ U^^^f-" 



2 V 2(2^ - k) - 2k - K 



In this case the self-interaction potential is asymptotically positive (see fig. 3). 

In the remaining cases, i. e., when either ^ = — ^ < (^ > 0), or k/2 < ^ < k, the 
resulting scalar field is an imaginary quantity so that, since we are interested only in real 
scalar field configurations, then the corresponding solutions will not be considered here. 

In this case, the scalar field solutions (7.13), and (7.15), both respect the 'mirror' sym- 
metry (5.1)-(5.5), but, unlike in the case 7.1, if is symmetric under x i— t- —x: v^^(x) = (p^{—x). 

7.3 General Case: ^ / 0, e / 

In this case we shall construct solutions that involve both the non-minimal coupling of 
the scalar field to gravity and the presence of the Gauss-Bonnet term. Despite the fact of the 
highly non-linear character of the relevant differential equations, we shall be able to obtain 
exact solutions that describe smooth braneworld configurations for this model. 
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Figure 3. V corresponding to solution (7.15). We set ^ = 10, ao = 1, 6 = 1 and k = 1/2. 



7.3.1 Particular solution: case u) = —1 

If one sets y = arcsinh x and = ^^^^^~'^\ then (4.4) can be easily solved for the 
particular case when oj = —1 or, equivalently, k, = 2^ since in this case we get a second order 
linear differential equation. Actually, by changing the variable x = sinhy, hence, cosh y = 
y/l + x"^, we get the following relevant equation for ip: 



2 ^ V «o 



((/3 + tanhy 99) — (J (99)^ — - sech^y (/?^ = - sech^y ( — 2~tanh^y— l). (7.16) 



By further performing the following substitution (f = ^^/(^ we get 

'ip + tanhy tp — ^ {1 — oj) sech^y il; = ^ {1 — u) sech^y [k^ tanh^ y — l) ip^^ , (7-17) 
2 s 

where k'^ = . By choosing the particular case u = — 1 we actually linearize this differential 
equation and get 

'i/' + tanhy ■0 — Ssech^y ip = — sech^y {k^ tanh^ y — l) . (7-18) 

This equation has the following real solution: 

2 14 - lO/c^ + e/s^sech^y ^ , , ^ ■ , , 
= 9? = — + Ci cosh {A) + C2 smh {A) , (7.19) 

where A = 2\/3 arctan (tanh |) , and Ci and C2 are arbitrary constants. By going back to 
the variable x we get the solution for the field 99: 



V9 = J ^ + ^^^^^^2) + ^1 (^) + ^2 sinh (A), (7.20) 



where now arctan 



\/x^+2-2Vx'^ + l 



This solution fulfills the symmetry requirements (5.1)-(5.5) when C2 = 0. Similarly to 
previous cases, V is asymptotically constant as it is shown in fig. 4. 
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V 




Figure 4. V associated to solution (7.20). We set Ci = 0, C2 = 0, ^ = 1/4, ao = 1, & = 1, e = 1/3 
and K = 1/2. 

7.3.2 Solutions for arbitrary uj 

Let us recall that if we perform the coordinate transformation x 1— )• r = arctan x (see Eq. 
(7.1)) for the case in which the Gauss-Bonnet term is non— trivial, then the master equation 
(4.4) adopts the form: 

^~'^\.-.2_3 2 3/ t2„- 2, 



which will be from now on the relevant differential equation to be solved. 

Let us now consider more general solutions in which the value of the parameter oj is 
arbitrary. When dealing with equation (7.21), we can perform the following transformation 
in order to reduce the order of this differential equation 

ip{T) = (^(T)/(r), (7.22) 

where /(r) is an arbitrary function of r. We further can redefine the scalar field if as follows 

ip^ = i^{T). (7.23) 

The relevant equation (7.21) transforms into 

Ifi^ + (/ + c^/' - ^ = -| (1 - fc' sin' r) , (7.24) 

With the aid of these two transformations, the second order differential equation (7.21) can 
be reduced to a system of first order differential equations, namely, to a non-linear first 
orderdifferential equation (general Riccati equation) for /(t) and a linear one for the new 
function iP{t): 

f + cof-^ = ^fH{r), (7.25) 

, , 3 (1 - fc^sin^r) 
i; + H{t)^ = ^ ^, (7.26) 
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where H{t) is a function of r. Solving this system is more easy than solving the original 
differential equation (7.21) if we make a suitable choice of the function H(t). 

By setting H{t) = 2ci//, with ci = const, in the general Riccati equation (7.25) we get 
the following solution for /: 



/ 



tan 



loo (r - To) 



(7.27) 



where 21 = — (2ci + 3) and tq is an arbitrary phase. We further substitute this solution into 
the linear equation (7.26) for •iP{t) and proceed to solve it. The general solution for this 
equation with arbitrary / possesses a lengthy expression and is given in terms of products 
of several hypergeometric and exponential functions combined with a sine at certain power. 
For the sake of simplicity we shall restrict to the case in which / = 1/w, getting the following 
solution: 



6uj [4 + 3a; - (2 + 3uj)k'^ sin^ (r - tq)] 
(4 + a;)(2 + a;) 



+ C2sin-(2+3-) (r-ro) 



(7.28) 



where C2 is an arbitrary constant. Once we have a solution for ip we can get back to the 
original field (p = according to (7.23). This solution will be subject to the relation (7.22) 
which implies that C2 = 0, leading to: 



6uj [4 + 3a; - (2 + 3u)k'^ sin^ (r - tq)] 
(4 + L^)(2 + a;) 



(7.29) 



Moreover, a further relation between k and oj completely determines the solution for ip. For 
instance, when fc^ = 2X3U) ' S^t a sine or cosine function depending on the phase constant 
Tq. It should be pointed out that both the original differential equation (7.21) and (7.26) are 
invariant under a suitable simultaneous rescaling of the field (p and the parameter ^. This 
fact that can be used to arrive at the following particular solutions for the field p: 



a) (f 




sin(r), 



k^ = l 



h) = ^/ — cos(t), 
5C 



k' 



2\uj\ 
2uj 



with CO < —5/2; 



1 H with 00 > -5/2. 

5 



(7.30) 
(7.31) 



In fig. 5 we show the profile of the V for the solution (7.31)^. Again, it is constant and 
negative at w — )• 00. There is another solution of the same type that is valid for the special 
value 00 = —11/8: 



c) ip = Bcos{2t) + C, 



(7.32) 



where the following constants 

B± = : 
C± = : 



154^ 



7 (2 - P) ± 2^4^4 + 49(1- A:2) 
33 k^ 



7 (2 - A;2) ± 2^4fc4 + 49(1 - k^) 



For solution (7.30) of the scalar field the qualitative behavior of V is similar. 
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Figure 5. V associated to solution (7.31). We set ^ = 1/3, ao = 1, 6 = 1, e = 1/5 and k = 1/2. 



must have the same sign (provided that the radicand is positive) since one can show that 
both parameters and k can be expressed as 



^ 11^2 + + 3C2 ' 

2 _ 28BC 
~ 11^2 + + 3C2 ' 

implying that the left hand sides of these equalities must be positive. 

One can look for more complex solutions than the ones displayed here. However, the 
relations that define the involved integration constants become more and more lengthy and 
difficult. 

From the above obtained exact solutions, only solutions (7.30) and (7.31) meet the 
symmetry conditions (5.1)-(5.5). 

8 Gravitational Fluctuations 

In order to study the localization properties of gravity, the analysis of gravitational 
fluctuations is needed. In 4D standard cosmology, the fluctuations of the geometry can be 
classified into scalar, vector and tensor modes with respect to the three-dimensional rotation 
group [57, 58]. This fact makes more feasible the study of the metric fluctuations because at 
the linear level the dynamical equations of the scalar, vector and tensor modes are decoupled. 
In our case, since the 4D geometry is Poincare invariant, the fluctuations of the metric may 
be also classified into scalar, vector and tensor modes with respect to the transformations of 
this symmetry group (see [59] for details). 

Let us consider the fluctuations of both the metric and the scalar field around the grav- 
itational background specified in (2.4) and the equations (2.5). In other words, the perturbed 
geometry has the following form 

dsl = {a^{w)r]AB + Has) dx^dx^, (8.1) 
while the fluctuation of the scalar fleld is 

H^P = V + X- 
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The functions Hab and x are the metric and the field fluctuations, respectively. The index 
p denotes perturbed quantities. 

By taking into account the 4D Poincare symmetry of our background metric (2.4), the 
fluctuations can be written as 



^AB - + ^AB + ^AB^ 



where 



Hf,=aHu,) l^no^'P + W , (8.2) 

+ (8.3) 

^fS-^H (8.4) 

The upper indices S, V and T denote the scalar, vector and tensor parts of the fluctuations, 
respectively. The tensor h^i, is transverse and traceless with respect to the 4D Minkowski 
metric ?7^jy, in other words 

= 0, duh"^ = 0. (8.5) 

Also, the vectors and Z)^ are divergence free 

d'^f^ = 0, di'D^ = 0. (8.6) 

The four remaining functions ip, E, C an are scalars with respect to 4D Poincare transfor- 
mations. 

The relations (8.2)-(8.6) tell us that, apparently, we only have 10 independent degrees 
of freedom of the metric fluctuations. Moreover, the covariance of our setup implies that the 
gravitational perturbation theory has some gauge degrees of freedom [57, 58]. In our case we 
can completely fix 5 of the above 10 degrees of freedom by making use of this gauge freedom. 
Let us choose what is known as the longitudinal gauge 

E = 0, C = 0, f^ = 0. (8.7) 

Thus, we finally have only 5 independent degrees of freedom. 

8.1 Tensor modes 

As in [60] the equation for the evolution of the tensor modes is 

- - -□"^Mv = 0, (8.8) 

where ^f^i, = s{w)hf^,^ with s{w) = a^q. In addition, r{w) = ^'^^^ ^ and denotes 

the d'Alembertian with respect to the metric ij. 

In order to study the mass spectrum of these fluctuations let us consider the following 
separation of variables ^^i, = '&{w)e^y{x). In terms of this new variable the equation (8.8) 
takes the form 

□"e- +m2e- =0, (8.9) 

- ^^dm + n?-^„, = 0, (8.10) 

yjs S 
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where e^(x) is a 4D tensor mode with mass m; on the other hand, i^miw) is the 5D profile 
of the field ^fj,i, and characterizes its localization properties. The equation (8.10) can be 
interpreted as a Sturm-Liouville eigenvalue problem with the associated norm: 

{^\^) = / -^'^dw. (8.11) 

J — oo ^ 

The zero mode e^^{x) is the massless 4D graviton. This mode is localized on the brane if 
(^ol^^o) = § ^0 finite. It is not difficult to show that ??o = \/s, then, the localization 
condition for the massless graviton is 

(^ol^o) = IZ a'qdw-Ae [a'^oo + 8 ^ ^oo ^ (8-12) 

As one can observe the above expression is quite similar to (3.5). Therefore, for the geometry 
described in (4.3), we have that if the 4D massless graviton is localized on the brane, then 
the 4D Planck mass is finite. 

Let us investigate the localization properties of the massless graviton in the backgrounds 
studied in section 7. The situation where = 0, e 7^ and k = 1 was analyzed in [61]. In 
this work it was shown that the zero mass graviton is localized on the brane. Furthermore, 
by substituting the expression (7.3) for A; / 1 into (8.12) we obtain the following result 

(^o|^o) = ¥(i + t)- (8-13) 

The above expression tells us that the 4D graviton is localized on the brane. 

In the following we shall consider the opposite case where ^ 7^ and e = 0. The 
normalization condition depends of the behavior of the integrand a?L{ip). With regard to 
the study of gravity localization, the backgrounds defined in (7.13) and (7.15) are similar 
because both have a regular and finite non-minimal coupling function L{ip). Therefore, the 
convergence of (8.12) is completely determined by the behavior when w — t- 00. For the 
warp factor (4.3) a'^ ~ l/|i(;p at infinity. This fact implies again that the massless tensorial 
mode is localized on the brane. 

Finally, let us study the general case (^ 7^ and e 7^ 0). Some exact solutions for this 
general situation are shown in (7.20) and (7.30)-(7.32). The norm for the zero mass mode 
takes the following form 

/oo 
{L-l)dw + {M^o)GB, 
'OO 

where {'&o\^o)gb is the norm written in (8.13). On the other hand, the function L associated 
to (7.20) and (7.30)-(7.32) is finite along the extra dimension, thus, like in the previous cases 
the normalization condition depends on the behavior of at infinity. From this fact we can 
conclude that the zero tensor mode is also normalized. 

As we mentioned above, in our model the graviton localized on the brane is synonymous 
of a finite 4D Planck mass. But, what about with the remaining coupling constants of (6.3)? 
The constants Ai and A2 are bounded due to the functional dependence of the warp factor 
with respect to w (4.3). Additionally, for all our solutions, \(p\, \(p'\ and \V\ are finite along 
the extra dimension and the self-interaction potential is asymptotically constant, thus, the 
behavior of the warp factor determines whether or not the effective 4D cosmological constant 
A4 is well-defined. The warp factor (4.3) has the following form at if —t- 00 

1 

a ~ -. — 7. 

\w\ 
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This fact ensures the convergence of the integral that defines A4 for all the solutions presented 
in section 7. In other words, the effective cosmological constant is finite. 



8.2 Vector modes 

In contrast with the tensor sector, the simultaneous presence of the Gauss-Bonnet term 
and the non-minimal coupling interaction makes much more difficult the study of vector 
modes. Thus, in this work we will study each case separately and will leave the general case 
for a future investigation. 

8.2.1 Vector modes with the Gauss Bonnet term only 

In this subsection we analyze the case where the non-minimal effects are negligible 
= and e ^ 0). Using (8.1), (8.3), (8.6) and substituting them in (2.2) we obtain the 
equations for the vector modes on the longitudinal gauge 

(Z)*^)' + (^^ + 3^^ D'^ = 0, (8.14) 
U'^B,, = 0, (8.15) 

where the first relation is a constrain and defines the profile of along the extra dimension. 
Furthermore, since there is no mass term in the above second equation, it shows that there 
is only a massless mode, the graviphoton, with no massive vector fiuctuations. 

For the case of vector modes, the issue of defining the norm is more indirect than for the 
tensor sector, therefore it is more appropriate to use the perturbed version of the action (2.1) 
up to second order with respect to the vector fiuctuations [54, 61]. This perturbed action can 
be written as follows 

,5(2)5^ = J d^xdw^ (r?"^5«P^5/3p^) , (8.16) 

where 2?^ = a^^'^^D^ is called canonical normal mode. The zero mode associated to (8.15) 
is localized on the brane if 5^'^^ Sy is finite. By making a suitable variable separation, the zero 
mass mode takes the form 



a''v''{x) = 0, 

where v^{x) is the 4D part of the vector sector of fluctuations. 
The norm of this mode reads 

(X,nP,.>^£5^. (8.17) 

When studying the localization properties of the zero mass tensor mode, we learned that it is 
necessary to have a convergent integral of a^q. Furthermore, the norm of the massless vector 
mode has an opposite behavior compared to the tensor mode one. Therefore, if one wishes 
the 4D graviton to be localized on the brane, the vector sector (represented by its zero mass 
mode) will necessarily not be localized on it. 
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8.2.2 Vector modes with the non minimal coupling only 

In this subsection we shall study the situation when e = and ^ 7^ 0. One way to obtain 
the mass spectrum of the vector fluctuations consists in passing from the action (2.1) (with 
e = 0) to the Einstein frame and then, perform the perturbation analysis. In order to do that 
we apply a conformal transformation as follows 

9ab = L^'^QAB- 
In the new frame the background action can be written as 



S^Sef = hk d''x^\^-R + \{V<jf - y(a)|, 



.18) 



where all quantities with an overline are associated to the metric , a is the new scalar field 
and V is its self-interaction potential. The relationship between the old and new variables is 



V 



Moreover, the general fluctuations on the Einstein frame are related with the old quantities 
as follows 

Hab = —^xL^'HAB + L^'^Hab. (8.20) 

The first term in the above expression does not contribute to the vector sector since it belongs 
to the scalar sector, then it is easy to obtain that = D^. Thus, the equations for the vector 
modes are 

(d'')' + mo'' = 0, 

UW^ = 0, 

where H. = a! /a and a = aL^/^ . 

Similarly to the previous case, let us define a new vector variable = a'^/^D^. There- 
fore, the norm of the massless mode takes the form 

dw dw 



The localization properties of this case are similar to those of the previous one. The localiza- 
tion of the massless tensor mode implies the delocalization of the vector sector of fluctuations. 

8.3 Scalar modes 

In the longitudinal gauge the scalar fluctuations of the geometry can be expressed as 



Furthermore, the matter provides an extra degree of freedom x to the scalar sector. Hence, 
it follows that we have three independent scalar degrees of freedom. One can obtain the 
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dynamical equations of the scalar fluctuations by perturbing the Einstein and Klein-Gordon 
equations with respect to the scalar sector. These equations adopt the following form 



4gV" + 4^/^' 
+4C' 



+ nq 



3 dip 



+ip'x' + m'xL^ + 



1 



0. 



?.21) 



mi - ^ [2%' + m^) + ^ 



+2C 



in' +n^) + 2nip' 



■ ■ • ■ 3 



ixL^r-n^ixL 



+ 2nxL^y 

= 0, 



^22) 



qC-2ij 



L 



xLw = 0. 



^ + 3q{^l,' +nC,) + {xL^)' -UxL^ + L'C = 0. 



(8.23) 
(8.24) 



x" + mx' - a^x + ¥^'[4V' + C] + 2C(/ + 3?^¥.') 

"I^a^x - ^|2(2?^' + 3'H')xL^^ + 

-3n^V + 4C(2?^' + sn^) + 4« + ^KV'' + cT 



0. 



(8.25) 



The expressions (8.23) and (8.24) are constraint equations. Then, there is only one indepen- 
dent scalar degree of freedom. As in the vector sector, in this case we will study separately 
the effects of the Gauss-Bonnet term and the non-minimal coupling on the scalar modes. 

8.3.1 Scalar modes with the Gauss Bonnet term only 

Although some aspects of this case were studied in [61], it is helpful to consider some 
of its details. The master equation of the system can be obtained by using the constraints 
(8.23), (8.24) and the dynamical equations (8.21), (8.22) and (8.25): 



3/2 ^, ^ . 

where $ = '^ ■^ and z = ° 'y^ . In order to study the mass spectrum of $ let us assume 
that 



(8.26) 



Thus, the equation (8.26) can be written as 



.27) 
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The associated norm for the above eigenvalue problem is 

($1$) = y (^1 + ±-^,^^dw. (8.28) 

The scalar massless mode is localized on the brane if its norm is finite. In other words 

($o|^o) = j + ^) '^odw < oo, (8.29) 

where $o = ^- Let us apply this general analysis to the background solution (7.3). It is not 
hard to show that when w — )• oo 



\w\'^ , when k = 1 

(8.30) 

\w\^, when k ^ 1 



Therefore, in both cases the massless scalar mode is not localized on the brane. 

8.3.2 Scalar modes with the non minimal coupling only 

Similarly to the vector modes case, the analysis of the scalar sector without the Gauss- 
Bonnet term is more easy in the Einstein frame. In this case the metric fluctuations can be 
expressed as follows 

flS=2#(Y°,)' (8-=") 

where ip = ip -\- and C = C ~ ^^ti'- other hand, the scalar field fluctuation in the 

Einstein frame x is related to x &s follows 



1 8 /L, ^ 2 



By considering the above expression and the equations (8.21)-(8.25) when e = 0, we obtain 
the master equation for the scalar sector of fluctuations 



= ) $ + m^$ = 0, (8.32) 



1^" 



z 

where $ = ^^-i-^ and z = " — . Thus, the norm of zero mode reads 

($o|$o) = / <^>Qdw = / -^dw. (8.33) 

The localization properties of the solutions (7.13) and (7.15) are similar in the sense that 
their massless mode has the same behavior at infinity, i.e. 

1 , ,5 

zjy ^ \w\ , when w ^ oo. 
The above result tells us that the zero mass scalar mode is not localized on the brane. 
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9 Conclusions 



We present a scenario with a thick braneworld model builded by a scalar field non- 
minimally coupled to the Einstein-Hilbert term. Furthermore, there is a Gauss-Bonnet term 
in the bulk. We obtain the whole 4D effective theory coming from the dimensional reduction 
of the set-up (2.1). We make emphasis in the analysis of relevant quantities such as, 4D 
Planck mass and the effective cosmological constant. As consequence of the non-minimal 
coupling between the scalar matter and gravity Mpi and A4 depend explicitly on the profile 
of the bulk scalar field ip. Therefore, in order for us to be able to rigorously interpret our 
physical results we need to construct exact solutions for our braneworld configuration. Despite 
the highly non-linear nature of the relevant differential master equation for the scalar field, 
we were able to obtain several exact particular solutions, contrary to previous studies that 
implement approximate solutions and/or particular cases that consider just one of the two 
effects. In order to elucidate the physical effects of the inclusion of the non-minimal coupling 
and of the Gauss-Bonnet term, we considered three cases in which we switched on/off the 
corresponding parameters. For all the constructed solutions, the effective Planck mass and the 
4D cosmological constant are finite. We also studied the whole set of gravitational fluctuations 
- classified into scalar, vector and tensor modes with respect to the 4D Poincare symmetry 
group - for our braneworld model. For all the considered backgrounds the massless zero mode 
of the tensor fluctuations is localized on the brane. This fact allows us to recover 4D gravity 
in our world. In contrast to this, the massless scalar and vector modes are delocalized and 
hence, cannot be observed on the brane. 

As a further step it will be interesting to study the stability of the field configurations 
within the the framework of the present braneworld model when both non-linear effects are 
switched on. Another interesting topic is to consider a non-minimal coupling of the bulk 
scalar field not only with the Einstein-Hilbert piece of the Lagrangian but, also, with the 
Gauss-Bonnet term. Investigation of these issues is left for future work. 
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